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Abstract 

Reproducing with elementary resources the correlations that arise 
when a quantum system is measured (quantum state simulation), al- 
lows one to get insight on the operational and computational power of 
quantum correlations. We propose a family of models that can simu- 
late von Neumann measurements in the x — y plane of the Bloch sphere 
on rt-partite GHZ states. 

1 Introduction 

Understanding the nonlocal correlations created upon measurement of 
some entangled quantum system is a problem which runs up against 
our common representation of the world, by the very definition of non- 
locality, i.e. violation of a Bell inequality |T]. Indeed, no explanation 
one would reasonably accept as possible, like agreement prior to mea- 
surement, or subluminal communication of inputs, seems to be used 
by nature in order to create these correlations (see the numerous ex- 
perimental violations of Bell inequalities ) . 

Still, some insight on the power of such correlations was gained 
when people came out with models able to reproduce them in terms 
of classical resources. For instance. Toner and Bacon 31 showed how 
to simulate von Neumann measurements on a singlet state with one 
bit of communication. Such a result puts an upper bound on the 
required amount of nonlocal resources needed for the reproduction of 
singlet correlations; it guarantees also that such correlations are not a 
stronger resource of nonlocality than 1 bit of classical communication. 

A different kind of resources that was also considered are the so- 
called nonlocal boxes [4 : these are simple nonlocal correlations which 
don't allow signaling. Successful simulation schemes using nonlocal 
boxes as unique nonlocal resources include the simulation of the singlet 
[S] and of partially entangled two-qubit states [B] . 



Concerning multipartite systems, communication models reproduc- 
ing Pauli measurements on n-partite GHZ or on graph states have also 
been proposed [7J [S] . For arbitrary possible measurement on the tri- 
partite GHZ state, previous studies suggested that its simulation with 
bounded communication might be impossible, taking as an example 
correlations corresponding to measurements of this state in the x — y 
plane of the Bloch sphere 9J. In this paper, we construct a model which 
reproduces precisely these equatorial correlations, and whose only non- 
local resources are Popescu-Rohrlich (PR) boxes ITp] and Millionaire 
boxes [llj . Note that even though our model doesn't give an upper 
bound on the worst-case communication cost, it can be translated into 
a communication model with finite expected communication cost, sim- 
ulating for instance the tripartite GHZ state with an average total of 
10 bits of communication between the parties (c.f. Appendix B). 

The paper is organized as follows: first, we recall the correlations 
of the GHZ state that we want to simulate. We then present a model 
for the 3-partite case, and generalize it to more parties. We discuss 
the construction and then conclude. 

2 GHZ correlations 

Consider the n-partite GHZ state 

|Gi/Z„) = ^(|00...0) + 111...1)). (2.1) 

Our goal is to reproduce the correlations which are obtained when von 
Neumann measurements are performed on this state, by using other 
non-local resources such as non-local boxes (possibly supplemented 
with shared randomness). 

For n = 2, the protocol presented in [5] for the singlet state allows 
one to reproduce the correlations for any measurement settings, using 
one PR box. Going to n > 3, we shall only consider measurements in 
the x — y plane (equatorial measurements) , which have the nice feature 
of producing trivial marginals: all correlation terms involving fewer 
than n parties vanish. Writing each party's measurement operator 
as: A ~ cos (f) a <Jx + sin (pa lyy , B — cos ajf -I- sin 0f, ay , . . . and 
denoting the binary result of each measurement bya,/3,...e{— 1,1}, 
the correlations we are interested in are given by 

(a) = (/?) = ... = {ap) = . . . = (2.2) 

for all sets of fewer than n parties, and 

{a(3...uj)=cos{(j)a + (t)b + ■■■ + (i)z) (2.3) 

for the full n-partite correlation term. In other words, outcomes appear 
to be random except when all of them are considered together, in which 
case their correlation takes a form reminiscent of the singlet state. 

It is worth mentioning that even though we restrict the set of possi- 
ble measurements on the GHZ states, the correlations we consider can 
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Figure 1: Simulation of \GHZ3) in a Svetlichny scenario: Alice and Bob 
form a group and can share their information with each other, while Charlie 
is separated from them. In this scenario 1 PR box allows one to reproduce 
the equatorial correlations. 



still exhibit full n-partite nonlocality. Indeed, the Svetlichny inequality 
for n parties can be maximally violated with settings in the x — y plane 
[T5J[T3]. This implies that in order to simulate these correlations, any 
model must truly involve all n parties together |14j . 

3 Simulation model for the 3-partite GHZ 
state 

Let us consider the above correlations for n = 3 parties, for which the 
outcomes of all parties need to be correlated according to (a/37) = 

COs((^a + 0h + 0c)- 

As a first step towards the simulation of these correlations, let us 
relax some of the constraints and allow two parties to cooperate in 
a Svetlichny-like scenario [T5] (see Figure [ij : for instance Alice and 
Bob would be allowed to communicate with each other, but not with 
Charlie who is kept isolated from them. In such a scenario, the three 
parties could create correlations of the desired form with one PR box 
by using the protocol of [S] to generate outputs a and 7 that have a 
cosine correlation of the form (aj) = cos{4'ab + '/'c), with a fictitious 
measurement angle (pab = 4'a + 4>b- By then setting either a = a, (3 = 
+1 or a = —a,/? = —1 (each with probability 1/2), and 7 = 7, they 
would recover the desired tripartite correlations (a/37) = cos(0q + 0b + 

0c). 

Of course, letting Alice and Bob share their inputs is not satisfac- 
tory yet, as this would require signaling between them. We shall now 
see that it is actually possible to re-separate them, while keeping the 
tripartite correlation term unchanged. 

In order to do so, let us recall that the model used above to create 
the bipartite cosine correlation with a PR box works by asking the 
parties (here, Alice-Bob together and Charlie) to input in the box 
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Figure 2: Setup for the simulation of \GHZ^) in the x — y plane. The 
Alice-Bob group was split by using two C-boxes and a second PR box. 



terms of the form [S] 

X = Sg(cOs(0Qb - .^i)) + Sg(cOs((/>afc - (^2)) 

z = sg(cos(0c - ^+)) + sg(cos(0c - 

where V'lj ¥^2: V+j ^re hidden variables shared by all parties (that we 
shall define later), and where the sign function is defined as sg(a;) — 
if a; > and sg(a:) = 1 if a; < 0. One can see that in the Svetlichny 
scenario, Alice and Bob don't really need to share their measurement 
angles, but only the terms sg(cos((/)a + '/>h — '^\))- Hopefully, there is 
a way for Alice and Bob to compute this function nonlocally by using 
a Millionaire box (M-box). For convenience, let us define the following 
nonlocal box: 

Cosine box. A bipartite Cosine box (C-box) is a non-local box that 
admits two angles (f)a, 4>b € [0, 27r[ as inputs and produces locally random 
binary outcomes a,b G {0, 1}, correlated according to 

a + 5 = sg(cos((/)a + 0b)). (3.2) 



We show in Appendix A that a bipartite C-box is equivalent to 
a M-box. C-boxes are exactly what we need for our problem, as the 
following result shows: 

Result 1. Equatorial von Neumann measurements on the tripartite 
GHZ state can be simulated with 2 C-boxes and 2 PR boxes. 

Proof. The simulation can be realized with the following model; we 
refer to Figure [2] for the distribution of the non-local boxes between 
the three parties Alice, Bob and Charlie, and for the numbering of their 
inputs (denoted Xi , yi and Zi for each party respectively) and outputs 
(denoted 0^,6^ and q). 

Let Bob and Charlie share two independent random vectors Ai, A2 
uniformly distributed on the sphere ^2. We define A± = Ai ± A2 and 
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refer to 0+, V'- for their phase angle in polar coordinates. Let 

the parties input the following variables into their boxes: 

Xy = 0a, ^2 = <i>a, = fli + 02 

yi = 06 + <<5l, y2=0h + <P2, y3 = &l + &2 (3.3) 

zi = Z2 = sg(cos((/)c - <^+)) + sg(cos((/)c - <^_)). 
The three parties then output: 

a = (—1)"^, with A = ai + 

/3=(-l)^, with 5 = 61 + 63 (3.4) 
7 = (-1)'^, with C = ci+C2+ sg(cos((/)c - 

The output of each party is the XOR of outputs received from non- 
local boxes shared with all other parties. Since these boxes are no- 
signaling, a single output of any nonlocal box is necessarily random. 
The only way as not to get a random average correlation is thus to con- 
sider all parties together, since missing one produces a random term. 
All correlations involving fewer than 3 parties thus average to zero. 

Concerning the 3-party correlations, we have: 

A + B + C ^ {ai + bi) + {as + ci) -I- (63 -I- C2) -I- sg(cos(0c - 

= (fli + 61) + X3Z1 + ?/3Z2 + Sg(cos(0c - ip+)) 

= (fli + 61) + (ai + 61 + 02 + 62)21 + sg(cos((/)c - (p+)) 
= (fli + h) + sg(cos(0c - </?+)) 

+ [Sg{c0s{(j)a + 0b + <Pl)) + Sg(cOs((/>a + 06 + ^2))] 

X [sg(cos(0c - <f+)) + Sg(cos(0c - f-))] 
= Sg{Vab ■ Ai) + sg(c • A+) 

+ Sg{Vab ■ Ai) + Sg{Vab ■ A2) Sg(c • A+) + Sg(c • A_ ) 

(3.5)" 

where we defined Vab — (cos(— 0q — 0b),sin(— 0q — 0b), 0) and with 
c— (cos 0c, sin 0c, 0) being Charlie's setting. Following the proof of [S] 
(see also [3J ) , we find that the average of this quantity over the values 
of the hidden variables Ai and A2 is: 

(A + i? + C) = ^ ' • = ^ ~ + '^^ + '^-^ (3.6) 

which leads, as requested, to 

(q;^7) = COs(0a + 0b + 0c)- (3.7) 

□ 

Coming back to the Svetlichny construction we see that it was in- 
deed possible to split the Alice-Bob group by allowing them to share 
two C-boxes. Concerning the PR box, it had to be split also, into two 
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new PR boxes, in order to recover the desired result: the computation 
made by the PR box in the Svethchny setup is now performed nonlo- 
cally, by the 2 PR boxes, using inputs distributed over the 3 parties. 

We restricted here to measurements in the x — y plane, but with 
a slight modification, Charlie can actually choose any measurement 
basis. Indeed a way to understand the appearance of the model for 
the singlet state (the bipartite cosine correlation), in the Svetlichny 
scenario, is to realize that the fictitious measurement angle (j)ab — 4'a-'r 
(f)h that Alice and Bob used above corresponds to the direction in which 
they would prepare a state for Charlie if they were to measure their 
part of the original GHZ state in their respective bases. In other words, 
in the quantum scenario, when Alice and Bob measure the GHZ state, 
they prepare one of the two state 

|z±) = -^(|0)±e-(^»+«|l)) (3.8) 

for Charlie. But a way for them to prepare one of these states if they 
share a singlet (or rather, a bipartite GHZ or | <&"*") state) with Charlie, 
is by measuring their part of the |$^) state along (p^b, which is what 
they effectively do in our model. So in fact they prepare a state \z±) for 
Charlie, which he can measure in the direction he wants (in particular, 
outside the x — y plane). The only modification in the model needed 
for that is that Charlie should use zi = Z2 = sg(c • A+) + sg(c • A_) and 
C — ci + C2+ sg(c • A+) to allow his measurement to point outside the 
X — y plane. 

We do not claim that the above model is optimal. It could be 
that strictly fewer nonlocal boxes are actually enough to reproduce the 
same correlations. It is nonetheless remarkable that truly tripartite 
correlations can be simulated with bipartite nonlocal resources only. 

It is also quite surprising that the model we presented here does 
not need more shared randomness than in the bipartite case. It might 
possibly be that a model that would use fewer nonlocal resources would 
require more shared randomness. 



4 Simulation model for the 4-partite GHZ 
state 



In the previous section we showed how to split (pab from equation (3.1 1 
into two phases (pm'Pb, in order to re-separate the group formed by Alice 
and Bob in the Svetlichny scenario. It is in fact similarly possible to 
split (/ic in order to have a total of 4 parties into play: 

Result 2. Equatorial von Neumann measurements on the j^-partite 
GHZ state can he simulated with 4 C-hoxes and 4 PR boxes. 

Proof. The simulation can be realized with the following model, analo- 
gous to the previous one; we now refer to Figure [3] for the distribution 
of the non-local boxes between the four parties Alice, Bob, Charlie 
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Figure 3: Setup for the simulation of \GHZ4) in the x — y plane. Notice that 
if any of the 4 parties is taken away (together with the non-local boxes it 
shares with the other parties), we recover a setup with 2 PR and 2 C-boxes, 
which corresponds to the simulation setup for \GHZ^), as in Figure^ 



and Dave, and for the numbering of their inputs {xi, yi, Zi and Wi) and 
outputs {ai,hi,Ci and di). 

Bob and Charlie still share two independent random vectors Ai, A2 
uniformly distributed on the sphere 5*2. With the same notations as 
before, let the four parties now input the following variables into their 
boxes: 

Xi=4>a, X2 = 4>a, X3 ^ X4 ^ ai + a2 

vi = <i)h + ipi, y2^(l)b + ^2, V3. = yi = bi + b2 

Zi= (f)c - Z2 = cj)c - Z3 = Z4 = Ci + C2 (4.1) 

Wi = W2 = cj)dj W3 = W4 = rfl + d2 

The parties should then output: 

a = (—1)"^, with ^ = fli + 03 + 04 
(3 = (-1)^, with B = 61 + 63 + 64 
7 = (-l)'^, with C = ci + C3 + C4 
6 = (-1)^, with L» = di + ^3 + ^4 



(4.2) 



For the same reason as in the tripartite case, all correlations of 
fewer than four parties vanish. For the 4-partite correlation term, the 
calculation of A+B+C+D is straightforward, following similar lines as 
in the tripartite case. It leads to a similar expression as in (3.5), except 
that c should now be replaced by Vcd — (cos((/)c + ^d), sin(0c + 0). 
This leads to the requested 4-partite correlation term: 

(a/37(5) =cos(0a + 06 + 0c + 0d). (4.3) 

□ 



Again, there is no claim of optimal! ty for the above model, but it is 
also remarkable that truly 4-partite correlations can still be simulated 
with bipartite nonlocal resources only, and no more shared randomness 
than for the bipartite case. 
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5 Going to more parties 



5.1 Possible extension of the model to any number 
of peirties 

In the last two sections, we showed how to construct models for the 

simulation of GHZ states involving n = 3,4 parties by splitting the n 
parties into two groups. Each group then had to calculate functions 
of the form sg(cos(S(^a; + ipx)) with for instance ^<pai = 4'a + <Pb, 
ip\ = ipi. Now, if we consider more parties, splitting them into two 
groups necessarily results in at least one of the groups having more 
than two parties. One could for instance have n — 1 parties on one side 
and 1 party on the other side. The sign function that each group has 
to calculate thus involves in general more than two phase angles. This 
motivates the definition of a generalization of the C-box to n parties: 

Multipartite Cosine box. An n-partite C-box is a non-local non- 
signaling box that admits n angles (fii £ [0, 27r[ as inputs and produces 
binary outcomes ai e {0, 1}, correlated according to 



The outcomes of the box are locally random. Also, all correlations 
involving fewer than n outputs vanish. 

Multipartite C-boxes allow one to generalize our model to the sim- 
ulation of multipartite GHZ state with any number of parties, by sep- 
arating the n parties into two groups, consisting of k parties on one 
side and n — k parties on the other: 

Result 3. Equatorial von Neumann measurements on n-partite GHZ 
states can be simulated with 2 k-partite C-boxes -h 2 {n — k)-partite 
C-boxes + k{n — k) PR boxes (for any < k < n). 

Sketch of the proof. Following the previous constructions, the group 
with k parties needs to calculate nonlocally two terms of the form 
sg(cos(^i -h . . . + 0fc + fx)), which can be done by using two fc-partite 
C-boxes, and the other group can similarly do its job with two (n — k)- 
partite C-boxes. As it was the case for the 4-partite case, each party 
from the first group also needs to share a PR box with each other 
party in the second group. We thus understand that by separating 
the n parties into these two groups, a total of 2 fc-partite C-boxes + 2 
(n — fc)-partite C-boxes -|- k{n — k) PR boxes is sufiicient to simulate 
the correlations of the n-partite GHZ state measured in the x — y plane. 
Interestingly again, no more shared randomness than for the bipartite 
case is required. 




(5.1) 



□ 
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5.2 A simpler model 



If we allow the parties to share nonlocal boxes involving more than two 
parties, then there is actually a simpler model which uses a single n- 
partite C-box to reproduce the equatorial GHZ correlations (as defined 



by (2.21 and (2.31): 



Result 4. Equatorial von Neumann measurements on n-partite GHZ 
states can be simulated with a single n-partite C-box. 

Proof. Consider indeed the following strategy: Alice generates a ran- 
dom variable (p\ G [— 7r/2,7r/2] according to the distribution p{ip\) — 
I cos{ipx)- She inputs (t>a + fx in the n-partite C-box, while all other 
n — 1 partners simply input their measurement angle. From the out- 
puts a,b, . . . of the box, each party can compute the final outputs 
a = (—1)°,/? = (—1)^,... All correlations between the outputs that 
involve fewer than n parties vanish, while for the n-partite correlation 
term, they get, as requested: 

(a/3 ... w) = / (-1) ^s('°<^^+^^+'t>^+-+<^^))p(^^^)d^x 

J -it/2 

_\\ ^"^^ COS ifx dipx - \ cos Lpx dfx if < < tt 

\-\ I-n/2'^'''' cos ipx difx + \ Jli%_j:^^ cos ipx dipx else 

= COs{(f>a <j)b + ■ . ■ + (pz) 

(5.2) 

□ 



Note that in the bipartite case, this model gives a new, simplified, 
way of simulating the equatorial correlations of the singlet state with 
a single Millionaire box. It is worth noting that it does not require 
any shared randomness. It uses however a strictly stronger non-local 
resource than the model with one PR box [5], since M-boxes cannot 
be simulated with one PR box (c.f. Appendix B). 

Compared to this last simple model, our previous construction al- 
lows one to reduce the multipartiteness of the nonlocal boxes used to 
simulate the same correlations. Finitely many nonlocal boxes involving 
no more than [n/2] parties are sufficient to reproduce n-partite equa- 
torial GHZ correlations. In particular, for n < 4, bipartite resources 
are sufficient. 

If one really wants to use only bipartite nonlocal boxes, we show in 
Appendix C that multipartite nonlocal boxes with continuous inputs, 
binary outputs, and only fully n-partite non- vanishing correlations, can 
always be simulated with bipartite boxes, as it is the case for boxes 
with a finite number of inputs |16j . However, the construction we use 
is quite special, as the boxes we need can have inputs or outputs that 
cannot be written as real numbers. 
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6 Conclusion 



We proposed models reproducing the correlations of the tripartite and 
4-partite GHZ states measured in the x — y plane, with a finite number 
of bipartite nonlocal boxes. Extending our results to n-partite GHZ 
states was possible after releasing the requirement that the nonlocal 
boxes had to be bipartite. 

We believe that our results give a new motivation for finding whether 
or not the GHZ correlations can also be simulated in a bounded com- 
munication scheme. Note that our models can be translated into finite 
expected communication schemes, since PR boxes can be replaced by 
1 bit of communication and M-boxes (bipartite C-boxes) by 4 bits in 
average, as we show in Appendix B. This gives a model with an average 
of 10 bits of communication between the parties. Note that this model 
with finite expected communication could also be recast as a detection 
loophole model. 

More generally, it would be interesting to know whether the sim- 
ulation of n-partite GHZ states can always be achieved with a finite 
amount of bipartite resources only (for instance a finite number of M- 
boxes). Considering also measurements outside of the x—y plane seems 
quite challenging because of the non-trivial marginals that appear, but 
would certainly be of interest too. 

Finally, it would be worth studying other multipartite quantum 
correlations. The W state for instance, seems to be a good candidate 
for this, when measurements are again restricted to the x — y plane, 
because of the simplicity of its correlations. Indeed they only consist 
of bipartite correlation terms of the form {a(3) = - cos(0a — 4>b), all 
other correlation terms being for any number of parties 
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Appendix A. The bipartite Cosine box is equivalent 
to a Millionaire box 

Here we show that in the bipartite case, the Cosine box is equivalent 
to a MiUionaire box, up to local operations on the inputs and out- 
puts. The general n-partite C-box can thus also somehow be seen as a 
generalization of a M-box to more parties. 

It is indeed clear that a bipartite C-box is equivalent to a "sine 
box", that would take two angles (pai 4'b € ["TI"; '''[ as inputs and would 
output two locally random bits a,b £ {0, 1} with correlations satisfying 

a + b^ sg(sin(0Q - (f>b)). (7.1) 

Now, if (j)a G [^TT, 0[, Alice can input 0a + G [0, 7r[ in the sine box 



instead of (pa, and flip her output so that (7.11 is still satisfied. This 
also holds for Bob; we can thus assume that (pa, (pb G [0, 7r[. 

In that case, sg(sin(0a — 4>b)) = sg(0a — <l>b)- The sine box thus 
compares the values of the two real numbers 4>a,(t>b', this is exactly what 
a M-box would do! 
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Figure 4: How to realize a bipartite Cosine box from a Millionaire box. 



More precisely, to construct a C-box from a M-box (that takes two 
inputs x,y e [0,1] and whose outputs satisfy a + b = sg(x-y)), Alice 
and Bob can input x = ^{(l)a mod tt) and y = ^(— ((^b+7r/2) mod tt). 
From the outputs a and b of the M-box, they can calculate a = a+ [^J 

and 6 = b + [— '^''"^^^^ J , which satisfy a + b — sg{cos{(t)a + (f>b)), as 
requested. This construction is illustrated on Figure |4] 

Reciprocally, the M-box can trivially be reproduced with a C-box, 
if Alice inputs <pa = x and Bob inputs — —y — it 12. 

Appendix B. Expected communication cost of a Mil- 
lionaire box 

We show in this Appendix that a Millionaire box cannot be simulated 
with finite communication. We propose however a scheme to simulate 
it with 2-way communication, unbounded in the worst case, but with 
a finite expected number of bits. 

Suppose first that Alice and Bob have a finite number 2^ of possible 
inputs (k bits). We show with a crossing sequence argument '18] that 
any communication scheme that can simulate the outputs of a M-box 
for this number of possible inputs necessarily uses at least k bits. This 
shows that in the limit of infinitely many inputs (ie, for the general 
M-box), unbounded communication is required. 

Indeed suppose that a scheme using ko < k bits of (possibly 2- 
way) communication can simulate the M-box with 2*^ possible inputs 
on each side. In particular it can simulate it when Alice and Bob use 
the same inputs x — y. But since the number of all possible identical 
inputs (2^^) is greater than the number of possible messages exchanged 
by Alice and Bob during the communication procedure (2*^"), there 
must be at least two different pairs of identical inputs xq = y^ and 
xi — yi (with xq 7^ xi) for which the communication pattern is the 
same. This communication pattern is then also the same if the inputs 
are x = xq, y — yi ot if x — xi, y — yo, because Alice and Bob will not 
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see any difference. So if the simulation of the M-box produces outputs 
saying that Xq < yi — Xi, it will also say that xi < i/q ~ xq, which 
contradicts the fact that xq xi. 

Thus simulation of a M-box with 2*^ possible inputs on each side 
necessarily needs at least k bits to be exchanged between the parties. 
So in the limit k ^ oo, the required amount of exchanged bits goes to 
infinity. 

Here is however a simple model that uses a finite average of 4 bits 
of communication (2 bits in each direction) to simulate a M-box. 

Let us write the two inputs x,y G [0, 1] of the M-box in basis 2, so 
that each digit is either a or a 1. Consider the following protocol, 
starting with n = I: 

1. Alice sends her n^^ digit to Bob. 

2. Bob compares the bit he received with his n^^ digit and answers 
if they are the same and 1 if they are different. 

3. If Alice receives a 0, she iterates n and goes back to step number 
1. If however she receives a 1, then they both know which one of 
them has the largest input number. Alice can output a predeter- 
mined random bit, and Bob a bit correctly correlated to Alice's, 
so as to reproduce the behavior of the M-box. 

The average number of communication cycles needed in this scheme 
depends on the probabilistic distribution of x and y. In particular if 
these distributions are independent and uniform on the interval [0, 1], 
the probability that the protocol stops at the n**^ step is 2~", and 
therefore the expected number of rounds is 

Since each round uses 2 bits of communication (one in each direction) , 
a total of 4 bits of communication is needed on average. 

Similar ideas can also be used to simulate n-partite C-boxes with 
finite expected communication. 

Appendix C. Simulation of n-partite correlations with 
bipartite nonlocal resources 

Consider an n-partite probability distribution for continuous inputs 
£ K and binary outputs — {0, 1}, which contains vanishing corre- 
lations for up to 71 — 1 parties. We show that it can be simulated with 
only bipartite nonlocal boxes. This can be seen as a generalization of 
|16j , in which a similar decomposition was constructed for distributions 
with finitely many inputs in terms of PR boxes. 

To show this result it is sufficient to concentrate on the total cor- 
relation term X^j'** = fi^i) involving all parties' outputs a^, because 
all other correlation terms can then be put to zero by letting all pairs 
of parties decide randomly to permute their output or not. 
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Consider thus the n-partite correlation term. We proceed by re- 
cursion: starting with the case n = 2, in which it is obvious that 
any bipartite no-signaling correlation can be produced by a bipartite 
nonlocal box satisfying 

ai + 0.2 = f{xi,X2). (7.3) 

Now let us suppose that we have a model which can reproduce 
any correlation term for rt — 1 parties. Any rt-partite correlation term 
can then be simulated the following way: for each value z that the n*^ 
party's input x„ can take, define the following function for the n—1 first 
parties: fz{xi, . . . ,x„_i) = /(xi, . . . ,a;„_i,z). Each of these functions 
can be implemented by the scheme reproducing the (n — l)-partite 
correlation functions. Now each of the n—1 parties can collect the 
outputs ai{z) it received for each possible value of z, and plug them 
into a special kind of bipartite nonlocal box it shares with the last 
party. This box takes as inputs a function f : [0, 1] {0, 1} on one 
side (i.e. a continuous number of binary inputs) and a real number 
X G [0, 1] on the other side, and produces binary outputs a and b, such 
that 

a + b = f{x) (7.4) 

If all parties input f = ai{z) into such a box they share with the 
last party, and this last party inputs x„ into all of these boxes, writing 
the outputs of each of these boxes and a^, we can set the last party's 
output to be a„ — to get a correlation term: 

n n — 1 n—1 n — 1 

^fli = ^ai + a„ ^(a, -l-a^) = ^a,(a;„) 

i—l i—1 i—1 i—1 ^ ' ' 

— ^ ("^l ; ■ ■ ■ 5 '^n— l) — y ('^l 5 ■ ■ • ; '^n) 

as required. 

Note that this construction needs continuously many nonlocal boxes. 



To avoid that, one could collect all the boxes (7.3 1 that calculate 
/z {xi , Xj ) for all z into a single one that would output all the values at 
the same time. Note however that such a box would actually output 
continuous outputs of cardinality H2 (i.e. binary functions defined on 



M). Note that the other boxes (7.4) also admit such inputs on one of 
their side. 



14 



